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Charge ordering accompanied by lattice distortion in quasi-two dimensional organic con- 
ductors 6>-(ET) 2 X (ET=BEDT-TTF) is studied by using an extended Hubbard model with 
Peierls-type electron-lattice couplings within the Hartree-Fock approximation. It is found 
that the horizontal-stripe charge-ordered state, which is experimentally observed in 0- 
(ET)2RbZn(SCN)4, is stabilized by the self-consistently determined lattice distortion. Fur- 
thermore, in the presence of the anisotropy in nearest-neighbor Coulomb interactions Vij, the 
horizontal charge order becomes more stable than any other charge patterns such as diagonal, 
vertical and 3-fold-type states. At finite temperatures, we compare the free energies of various 
charge-ordered states and find a first-order transition from a metallic state with 3-fold charge 
order to the insulating state with the horizontal charge order. The role of lattice degrees of 
freedom in the realization of the horizontal charge order and the relevance to experiments on 
#-(ET) 2 X are discussed. 
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Quasi- two-dimensional molecular conductors (ET) 2 X 
(ET=BEDT-TTF) show a variety of electronic phases 
at low temperatures. 1 ' 2 Among them, charge order (CO) 
phenomena are one of the main subjects and have been 
intensively studied recently. (ET) 2 X is a member of the 
so-called 2:1 salts, which consists of alternating layers 
of anionic X - and cationic ET +1 / 2 whose 7r-band is 
3/4-filled. The variety of their physical properties origi- 
nates from the spatial arrangements of ET molecules and 
strong Coulomb interaction among tt electrons. 

The experimental observations of CO are made 
in compounds such as #-(ET) 2 RbZn(SCN) 4 3,4 and 
a-(ET) 2 I 3 . 5,6 <9-(ET) 2 RbZn(SCN) 4 shows a metal- 
insulator transition at T = 200 K and a spin gap be- 
havior at low temperatures. 7 The transition is of first 
order accompanied by lattice distortion. CO formation 
below T c has been directly observed in NMR experi- 
ments. 3,4 Several experiments 8 11 such as Raman scat- 
tering 9 and X-ray scattering 10, 11 measurements indicate 
that the horizontal-type CO is formed in this compound. 

CO phenomena are considered to be a consequence of 
strong correlation effects among electrons, especially due 
to the long-range component of the Coulomb interaction. 
So far, many theoretical investigations on CO have been 
carried out from this point of view. 12 19 For example, Seo 
studied the extended Hubbard model, that includes both 
on-site (U) and intersite (V) Coulomb interactions using 
the Hartree approximation and discussed the stability 
of various stripe CO patterns in (ET) 2 X. 12 For #-type 
salts, the possibility of CO with long periodicity has been 
considered within the Hartree approximation. 16 Actu- 
ally, X-ray experiments on #-(ET) 2 RbZn(SCN)4 indicate 
a short-range CO with long periodicity in the metallic 
phase which is different from the horizontal stripe state 
at low temperatures. 10, 11 A similar charge fluctuation is 

* E-mail address : yasuhiro@ims .ac.jp 



observed in #-(ET) 2 CsZn(SCN)4, which shows coexisting 
charge modulations with different wave vectors without 
long-range order. 20,21 

On the other hand, a coupling between electron and 
lattice degrees of freedom also seems to have an im- 
portant role. In fact, the CO transition is accompanied 
by a structural distortion in many #-type materials in- 
cluding #-(ET) 2 RbZn(SCN)4. Moreover, a particular role 
of structural modification at the transition is suggested 
by a recent observation of photoinduced melting of CO 
in 6>-(ET) 2 RbZn(SCN) 4 and a-(ET) 2 I 3 . 22 Several theo- 
retical studies 12,13,16 indicate that lattice effects indeed 
stabilize the horizontal CO in #-(ET) 2 X, although any 
electron-lattice coupling which causes structural change 
is not explicitly included in the calculations. Thus, it is 
important to investigate not only the role of electron- 
electron interactions but also lattice effects on CO. 

In this paper, we study the CO transition and lat- 
tice distortion in #-(ET) 2 RbZn(SCN)4 by using the 
extended Hubbard model with Peierls-type electron- 
lattice couplings within the Hartree-Fock approxima- 
tion. Figures 1(a) and 1(b) show the structures of 6- 
(ET) 2 RbZn(SCN)4 in the metallic and insulating phases, 
which is called #-type and #d-type, respectively. At T > 
T c , the unit cell contains two molecules and two kinds of 
transfer integrals, t c and t p . On the other hand, six trans- 
fer integrals exist in the unit cell with a doubled c-axis at 
T < T c . Since the displacements of ET molecules and the 
resulting change in transfer integrals are rather compli- 
cated, 10 here we study the effects of electron-lattice cou- 
plings which cause the modulations of transfer integrals 
that are experimentally observed [Fig. 1(b)], and do not 
consider any other electron-lattice couplings. This leads 
to three kinds of interactions between electrons and the 
lattice degrees of freedom: transfer integrals modulated 
by c- and a-axis molecular translations and rotation, as 
deduced from the results of the X-ray experiment. 10 For 
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Fig. 1. Schematic representation of (a) #-type and (b) O^-type 
structures. The gray solid lines indicate the unit cell. The values 
of transfer integrals in (b) from the extended Hiickel method are 
as follows, t c i = 1.5, t C 2 = 5.2, t v \ = 16.9, t P 2 = —6.5, t P 3 = 2.2, 
and t p 4 = -12.3(10- 2 eV). 



simplicity, these electron-lattice couplings are assumed 
to be independent of each other. First, the c-axis trans- 
lation alternates t c and gives t& and t C 2 in Fig. 1(b). 
This is indeed expected since the length of the t c ± bond 
increases while that of the t C 2 bond decreases through 
the CO transition. 10,23,24 On the other hand, we observe 
that the length of the t v \ bond decreases while that of 
the t p s bond increases, 10,24 from which the modulations 
of t p i and t p s can be regarded as due to the a-axis trans- 
lation. However, a similar consideration does not hold for 
the changes of t P 2 and t p ±. In fact, the experimental esti- 
mation of transfer integrals indicates that rotational de- 
grees of freedom are important. It shows that the depen- 
dences of the transfer integrals on relative angles (called 
elevation angles 10 ) of ET molecules are large and allow 
the transfer integrals \t P 2\ (\t P 4\) on the horizontally con- 
nected bonds to uniformly decrease (increase), as can be 
seen from Fig. 1(b). In the actual compound, this type of 
modulation seems to be important since the horizontal 
CO is formed by the t P 4 chains with hole-rich molecules 
and the t P 2 chains with hole-poor molecules. Therefore, 
in the present study, we simply introduce such rotational 
degrees of freedom in order to take account of the exper- 
imentally observed modulations of t p 2 and t p ±, which are 
difficult to understand from a molecular translation. 
Our Hamiltonian is then written as 



(1) 



i (ij) (ij) 

where (ij) represents the summation over pairs of neigh- 



boring sites, c\ a (ci a ) denotes the creation (annihilation) 
operator for an electron with spin a at the ith site, 
= c lcr c i<j, and rii = + n^. The transfer integral 
tij means t c or t p in the #-type structure. The electron- 
lattice coupling constant, the lattice translational or ro- 
tational displacement and the elastic constant are de- 
noted by aij, Uij, and Kij, respectively. The electron 
density is 3/4-filled and we consider nearest-neighbor 
Coulomb interactions V c for the vertical direction and 
V p for the diagonal direction as shown in Fig. 1(a). For 
the lattice degrees of freedom, we further introduce new 
variables as jjij = ctijUij and Sij = afj/Kij, where 
Si j is written as s c , s a and for c-axis translation, a- 
axis translation and rotation, respectively, as discussed 
above. Similarly, we can rewrite yij by using the sub- 
scripts c, a and </>, and as a result the transfer integrals 



(a) 3fold 



(b) diagonal 





Fig. 2. Order parameters considered in the Hartree-Fock approx- 
imation. The hole-rich and -poor sites are represented by the 
solid and open ellipses, respectively. 



in the distorted structure are given by 

t c i = t c + y c , t C 2 = t c — Vc , 

t P i = t p +y a , t P 2 = t p — y { f) , (2) 

tp3 — tp — y a •> tp^ = t p + y§ . 

Note that the signs in eq. (2) are chosen so that the 
resulting transfer integrals correspond with the experi- 
mental ones if y\ > for I = c, a and (j)? b 
We apply the Hartree-Fock approximation, 



(3) 



to eq. (1) and diagonalize the obtained Hamiltonian in 
fc-space by assuming the unit cell of each mean-field or- 
der parameter. We use four types of CO order param- 
eters with respect to charge degrees of freedom, which 
are shown in Fig. 2. As for the spin degrees of freedom, 
we consider three spin configurations in each stripe-type 
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CO which are identical to those of ref. 12. For the 3- 
fold CO, spin alternation between the hole-rich and -poor 
sites is considered. The ground-state energy is calculated 
by solving the mean-field equation self-consistently with 
the lattice displacements, which are determined by the 



Hellmann-Feynman theorem ^f^y = 0, where y\ means 
y c , y a or y^. The energy per site is given by 



Zkcr 



5> 

(ij) 



[ni)(nj 



(4) 



(ij) 



where Z, Ei\^ a and rip are the band index, the energy 
eigenvalue of the mean-field Hamiltonian and the Fermi 
distribution function, respectively. TV is the total number 
of sites. In the following, we set t p = 0.1 eV, t c = —0.04 
eV, and U = 0.7 eV. The ratio V c /U is fixed at 0.35 and 
the anisotropy in nearest-neighbor Coulomb interactions 
Vp/V c is treated as a parameter. 

The ground-state energies of various CO patterns per 
site as a function ofV p /V c are compared in Fig. 3, where 
the energy of the 3-fold CO is set at zero. We have shown 
only the lowest-energy state of each CO pattern with dif- 
ferent spin configurations. In the absence of an electron- 
lattice coupling, the 3-fold CO with a ferrimagnetic spin 
configuration is the most favorable in the nearly isotropic 
region, i.e., V p /V c ~ 1. On the other hand, the diagonal 
CO whose spin configuration is antiferromagnetic along 
the stripe and between stripes on the c-axis is stable 
when V p /Vc is small. These features are consistent with 
the previous study. 16 For the horizontal CO, we plotted 
the energy of the state which is antiferromagnetic along 
the stripe and ferromagnetic between stripes on the c- 
axis. Note that the state which is antiferromagnetic on 
the c-axis has a close energy and is nearly degenerate 
with the above state. As can be seen from Fig. 3, there is 
no region where the horizontal CO has the lowest energy 
in the absence of an electron-lattice coupling. 

However, in the presence of the electron-lattice cou- 
plings, the horizontal CO becomes more stable owing 
to the lattice distortions. The values of the electron- 
lattice couplings are chosen at s c = 0.08, s a = 0.17 and 
= 0.11 to obtain realistic values of lattice displace- 
ments. The horizontal CO has hole-rich sites on the t P 4 
chains, which is consistent with the experiments. The en- 
ergy gain mainly comes from the difference between t P 2 
and t P 4. This is reasonable since the horizontal CO can 
be stabilized by the exchange coupling between neighbor- 
ing spins on the stripes. Although the diagonal CO is not 
affected by any electron-lattice coupling, the vertical and 
3-fold COs also have energy gain from the lattice modu- 
lation. The energy of the vertical CO is lowered by s a . On 
the other hand, that of the 3-fold CO is lowered by s c and 
Sep. In this state, there is a weak horizontal charge modu- 
lation caused by the lattice distortion in the background 
of the 3-fold CO. Note that the 3-fold CO is metallic even 
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Fig. 3. (Color online) Relative energies as a function of anisotropy 
Vp/Vc, where the energy of the 3-fold state is chosen to be zero, 
dg, hr and vt are abbreviations of diagonal, horizontal and ver- 
tical COs, respectively. 
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Fig. 4. (Color online) (a) Modulations of transfer integrals and (b) 
hole density at each site in the horizontal CO. In (b), the hole 
densities in the absence of an electron-lattice coupling are also 
shown. 



if the transfer integrals are modulated, whereas the hor- 
izontal CO is insulating. As a result, the horizontal CO 
with lattice distortion becomes stable for V p /V c < 0.75, 
while the 3-fold CO is favorable for V p /V c > 0.75, as 
shown in Fig. 3. 

In Fig. 4, we show the modulations of the transfer in- 
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tegrals and the hole density at each site in the case of the 
horizontal CO. Although y c and distortions do not de- 
pend so much on V p /V c , the y a distortion increases with 
V p /V c . It becomes largest for V p /V c > 0.7. In fact, the 
former two electron-lattice couplings favor the horizontal 
CO while the latter tends to decrease the order parame- 
ter as seen from Fig. 4(b), although the energy is lowered. 
Experimentally, the difference between t p i and t p s is the 
largest while that between t c i and t C 2 is the smallest. Our 
result seems to be consistent with the experimental one 
for V p /V c ~ 0.7. The detailed role of each electron-lattice 
coupling on the horizontal CO is discussed elsewhere. 26 
Next, we consider the stability of these COs at finite 
temperatures by calculating the free energy within the 
Hartree-Fock approximation. The free energy per site is 
written as 

F = Jj ^ Ntot ln(1 + e M-P(EiK* - A*)}) 

' ;k<T 

- U Y,( n ^( n n) -J2Vij(ni)(nj) 

+ E^< c l^)(sW) + E^) > 

(ij)<r (ij) %J 

where /i, N to t and (3 are the chemical potential, the total 
number of electrons and the inverse temperature, respec- 
tively. The phase diagram on the (T, V p /V c ) plane, which 
is obtained by comparing the free energies of different CO 
patterns, is shown in Fig. 5. The values of the electron- 
lattice couplings are the same as those used for T = 0. 
For Vp/V c ~ 1, the 3-fold state with lattice distortion 




v p /v c 

Fig. 5. Phase diagram on (T,V P /V C ) plane in the presence of 
electron-lattice coupling. 

has the lowest free energy for a wide temperature range. 
On the other hand, the horizontal CO is the most stable 
for Vp/Vc < 0.7. There is a first-order metal-insulator 
transition from the 3-fold CO to the horizontal CO near 
V p /Vc — 0.7. For the parameters we used, the normal 
metallic state with a uniform charge density always has 
a higher free energy than the 3-fold CO because of the 
large Coulomb interactions. We note that if we choose 
smaller values of U and Vij, the diagonal CO, which is 



undistorted even with electron-lattice couplings, is more 
stable than the horizontal CO. Therefore, large Coulomb 
interactions seem to be important to stabilize the hori- 
zontal CO with realistic values of the lattice distortions. 
In fact, the energy gain of the horizontal CO due to the 
lattice distortions can be understood by the perturba- 
tional calculations from the strong couping limit, i.e., 

Us = o. 27 

Finally, we discuss the relevance of the results to the 
experiments and relations to other theoretical studies. 
The stabilization of the horizontal CO due to the lat- 
tice distortion is consistent with the experiments on 0- 
(ET) 2 RbZn(SCN) 4 . Since the horizontal CO does not 
become the ground state without electron-lattice cou- 
pling, the effects of the lattice distortion are considered 
to be crucial in realizing the horizontal CO in the present 
model. This result is also qualitatively consistent with 
the recent exact-diagonalization study 27 for eq. (1) on 
small clusters. Moreover, the first-order metal-insulator 
transition at a finite temperature can be related to the 
experimental results of this compound, although the ob- 
tained wave vector of the charge modulation at high tem- 
peratures is different from that of the experiments in 
the metallic phase. It has recently been pointed out that 
longer range than the nearest-neighbor Coulomb interac- 
tions can reproduce the experimental observation. 18 As 
for the spin degrees of freedom, both the 3-fold and hor- 
izontal COs in our Hartree-Fock calculation have spin 
orders which have not been observed in the experiments. 
It is considered that the effect of quantum fluctuation is 
necessary in discussing the behavior of the spin degrees 
of freedom. 16 

The previous estimations of the intersite Coulomb in- 
teractions V p and V c show that these values are compara- 
ble, Vp/V c ~ l, 14 where the 3-fold CO is the most stable 
in our calculation. A variational Monte Carlo study 17 
in the absence of an electron-lattice coupling also shows 
that the 3-fold CO is stable for V p /V c ~ 1. According to 
the recent exact-diagonalization study, 27 the horizontal 
CO with lattice distortion becomes more stable even at 
Vp/Vc °° 1 if we take account of quantum fluctuations 
that are neglected in the Hartree-Fock approximation. 
An exact-diagonalization study 13 also indicates that the 
Holstein-type electron-lattice coupling stabilizes the hor- 
izontal CO. 

At the nearly isotropic region V p /V c ~ 1, we find that 
the 3-fold state with a coexisting weak horizontal charge 
modulation is stable. This result can be related to the 
X-ray experiments on <9-(ET) 2 CsZn(SCN) 4 , 20 ' 21 which 
shows two types of COs coexisting as short-range fluc- 
tuations. Although the present Hartree-Fock calculation 
gives a long-range CO, it is natural to expect that the ef- 
fect of fluctuations can destroy the long-range order and 
results in a state such that two types of COs coexist as 
short-range fluctuations. 

In summary, we investigated the effects of Peierls-type 
electron-lattice couplings on the CO in #-(ET)2X by us- 
ing the extended Hubbard model through the Hartree- 
Fock approximation. We found that the horizontal stripe 
CO which is observed in the experiments is stabilized by 
the lattice distortion. Moreover, at finite temperatures, 
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there is a first-order metal-insulator transition in the 
presence of the anisotropy in Vij , which can be related to 
the phase transition in #-(ET)2RbZn(SCN)4. These re- 
sults show that the lattice effect plays an important role 
on the CO phenomena in 0-(ET)2X. 
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